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ABSTRACT 



Aims. We test the reliability of the observed and calculated spectral irradiance variations between 200 and 1600 nm over a time 
span of three solar rotations in 2004. 

Methods. We compare our model calculations to spectral irradiance observations taken with SORCE/SIM, SoHO/VIRGO and 
UARS/SUSIM. The calculations assume LTE and are based on the SATIRE (Spectral And Total Irradiance REconstruction) 
model. We analyse the variability as a function of wavelength and present time series in a number of selected wavelength 
regions covering the UV to the NIR. We also show the facular and spot contributions to the total calculated variability. 
Results. In most wavelength regions, the variability agrees well between all sets of observations and the model calculations. The 
model does particularly well between 400 and 1300 nm, but fails below 220 nm as well as for some of the strong NUV lines. 
Our calculations clearly show the shift from faculae-dominated variability in the NUV to spot-dominated variability above 
approximately 400 nm. We also discuss some of the remaining problems, such as the low sensitivity of SUSIM and SORCE for 
wavelengths between approximately 310 and 350 nm, where currently the model calculations still provide the best estimates of 
solar variability. 
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1. Introduction 

The solar irradiance, or the solar flux received at the top 
of the Earth's atmosphere, is known to vary over a large 
number of time scales, ranging from minutes to months 
and decades. The changes in the total solar output have 
been measured since 1978 (|Willson fc Hudsonlll988l ) and 
different c omposites of the measu r ements have been pre- 
sented by iFrohlich k LeanI (llQQSl) : IWillson fc Mordvinov 
(|2003l ) and iDewitte et al.1 (|2004l ). While the short-term 
(minutes to hour) variability is mainly due to solar os- 
cillations and granulation, the daily to decadal variability 
is attributed to the changes in the surface magnetic field 
combined with the solar rotation that trans ports solar ac- 
tive r egions into and out of view. Indeed, iKrivova et al.l 
(|2003h found that more than 90% of the solar variability 
between 1996 and 2002 could be explained by changes in 
the solar surface f ield. S imilar conclusions were reached 
by IWenzler et al. ( 20061 ) who reconstructed solar irradi- 



ance from Kitt Peak magnetograms covering the last 3 
solar cycles. 

Solar variability is a strong function of wavelength: 
while solar output is small in the UV, the relative vari- 
ability is more than one order of magnitude larger in the 
UV than in the visible. Until very recently, the spectral 
dependence of the solar variability had mainly been de- 
termined in the UV, in particular by the measurements 
taken by the SUS IM and SOLSTICE instruments onboard 
UARS (see, e.g.. iFlovd et al.l l2003al ) . Information in the 



visible was restricted to the three colour channels of the 



SPM instrument of SOHO/VIRGO (jFrohlich et al.lll995D . 
though degradation hampered the use of these data be- 
yond timescales of the order of a few monthfl 

The variability at most other wavelengths had to be in- 
fer red us i ng a v ariety of approaches, such as e.g., pioneered 
bv iLean ( 1989f) who produced the first estimate of solar- 
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|f tp: //f tp .pmodwrc . ch/pub/data/irradiance/virgo/SSI/spm_level2_d_17C 
see also.ftp : //ftp. pmodwrc. ch/pub/Claus/SDRCE_Sep2006/SSI_Poster.pdi , 
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cycle variability over a large wa velength rang e . An alter- 
native approach was followed bv lUnruh etlll (|l999l) who 
used facular and spot model atmospheres t o calculate the 
flux changes due to magn etic features. Fhg ge et al. ( 2000? ) 



and iKrivova et all ( 2003h used solar surface images and 
magnetograms to calculate the variability on time scales 
ranging from days to years. Here we built on this approach 
and present comparisons between modelled and measured 
spectral irradiances during three months in 2004. 

Thanks to missions such as SORCE and SCIAMACHY 
the observational outlook has now become much better 
and we have, for the first time, variability observations 



that span from the UV to the near IR ([Harder et al 



2005bl : iRottman etHI l2005t ISkupin et al]l2005l ). In the 

following we consider SORCE data only. First com- 
parisons between SORCE measurements and models 



have been presen ted by, e.g., iFontenla et al.l (|2004l ) and 
Lean et all (|2005[ ). 

All data presented here have been recorded between 21 
April and 1 August 2004. During this time the Sun was 
in a relatively quiet phase, especially in May when only a 
very small spot group appeared on the solar disk. A new 
and larger active region emerged over the next month, 
resulting in a depression of just over 1 permillc in total 
solar irradiance (TSI) in July. 

In the next section we briefly describe our irradiance 
modelling approach. We then discuss the data analysis for 
the different instruments (Sec. [3]). In Sec. IH we compare 
the relative irradiance changes derived from the models 
with a number of different data sets spanning a wavelength 
range from 200 to 1600 nm. In particular, we compare our 
model to data from SORCE/SIM, UARS/SUSIM, and 
SoHO/VIRGO. We conclude this section by presenting 
observed and modelled timeseries in a number of selected 
wavelength bands. A discussion of the results and conclu- 
sions are presented in Sec. [5l 



2. Irradiance reconstructions 

Here we restrict ourselves to a brief description of our ap- 
proach to model solar irradiance (see Fligge et al. 2000, 
Krivova et al. 2003 for a more detailed discussion). 
Essentially, we calculate the solar irradiance (or flux) by 
integrating over the (pixellated) solar surface, account- 
ing for the presence of dark (sunspots) and bright (facu- 
lae and network) surface magnetic features. The location 
of sunspots is obtained from MDI continuum images, at- 
tributing penumbra and umbra to those pixels with con- 
trasts of less than 0.9 and 0.6, respectively. Faculae and 
the network are identified by their excess magnetic flux 
on MDI magnetograms. As faculae are very small-scale 
features and typically do not fill an entire MDI pixel, we 
adopt a filling-factor approach, scaling the facular filling 
factor (linearly) with the magnetic field strength measured 
from the magnetograms. The identification of the mag- 
netic f eatu res is described more e xtensively in lFligge et al, 
(|2000[) and IKrivova et all (|2003l ). 



The model has a single free parameter, Bsat, which 
takes into account the saturation of brightness in re- 
gions with higher concentrat ion of magnetic elements 
(e.g., Solanki fc Stcnflo 1981 ISolanki fc Briglie-^ Il992 : 



Ortiz et^n 2002). i?sat denotes the field strength be- 



low which the facular contrast is proportional to the 
magnetogram signal, while it is independent (saturated) 
above that. From a fit to the VIRGO TSI time series 
Krivova et al. I (|2003h obtained a value of 280 G for i?sat, 



which is used here unchanged. 

Once each pixel on the solar surface has been identified 
as either (part) facula, quiet Sun, umbra or penumbra, we 
can attribute a corresponding emergent intensity to it and 
then proceed to carry out the disk integration. Note that 
the emergent intensities have to be known as a function of 
limb angle for each of the components present on the so- 
lar surface. The wavelength resolution and available range 
for the final spectral irradiance is determined by the wave- 
length resolution and range of the limb-dependent emer- 
gent intensities. 

We calculate the intensities from th e SATIRE set of 
model atmospheres ( Unruh et al. using Kurucz' 



ATLAS9 program (lKuruczlll993ir The model atmosphere 
for the faculae and network was derived using FAL P 
(£ontciila^ct_al. 1999) as a starting point, while the quiet- 
sun is Kurucz' standard solar atmosphere and the sunspot 
umbra and penumbra model s are ste l lar m odels at 4500 K 
and 5150 K also taken from Kurucz As the model 

atmospheres and intensities are derived under the assump- 
tion of LTE, we expect our irradiances to b ecome unreli- 
able below approximately 300 nm (see, e.g.. lUnruh et al 
1999HKrivova et alllioO 



For the comparisons presented in this paper, the avail- 
ability of MDI imagefH and groups of five consecutive mag- 
netograms (which were averaged to reduce the noise) was 
reasonably good and we were mostly able to calculate irra- 
diances on a 12-hourly interval. There are, however, some 
data gaps, most noticeably at the end of June with only 
three sets of images with poorer quality between 2004 June 
22 and June 30. 



3. Solar irradiance observations: May to July 2004 

The main instruments used for the comparisons are the 
spectral and total irradiance monitors from SORCE, SIM 
and TIM, respectively. These data are complemented by 
contemporaneous observations from UARS/SUSIM and 
VIRGO/SPM. In this section, we briefly describe the in- 
struments used and discuss the data analysis. 

SORCE (Solar Radiation and Climate Experiment) 
was launched in January 2003 and started science oper- 
ations in March of that year. It is the first satellite to 
provide reliable daily measurements of the spectral irra- 
diance variability for wavelengths longer than 400 nm. It 



^ All MDI images were obtained from the MDI homepage at 
|http: //soi ■ Stanford ■ edu/[ 
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ca rries four instruments , all of which have been described 
in 



Rottman et alj (j2005[ ) and sources referenced therein. 



0.008 



3.1. SORCE/SIM 

SIM primarily measures spectral irradiance between 300 
and 2400 nm with an additional channel to cover the 200 
to 300 nm wavelength region. We consider data taken 
with three of its five detectors, namely the UV detector 
(200-308 nm), the VISl detector (VISl: 310-1000 nm), 
and the IR detector (994 — 1655 nm). These will be dis- 
cussed briefly in the following sections. We discard the 
data from the second visible-light detector as it suffers 
from both temperature and radiation-induced variability 
that cannot be fully removed. We were unable to use 
the longer- wavelength data recorded by SIM/ESR as they 
were too noisy over the time span considered here. For 
more information on the design and calibration of SIM we 
refer to Hfu-dejr et al. ( 2005a , b). 

The results presented here are based on Version 10 
of the SIM data reduction. The availability of SORCE 
data in the time considered here is reasonably good with 
only some data gaps and correction problems during two 
weeks around the end of June and beginning of July. 
The SOHO/MDI suffered from poorer imaging data dur- 
ing these two weeks as well, making comparisons at these 
times more difficult. 

All three detectors provide measurements of the solar 
irradiance as a function of w avelength on approxim ately 
12-h intervals. As discussed in [Harder et all (|2005bh . SIM 
typically has 6 samples per resolution element, yielding 
an (un- aliased) oversampling by about a factor of 2. In 
order to compare the data to the model calculations, we 
characterise them in two ways. Firstly, we consider the 
variability, treating each wavelength bin as an indepen- 
dent time series. As a measure of the variability, we adopt 
the standard deviation, calculated according to 



(t(A,) 



in -I) 



(1) 



where f{Xi) is the mean flux at wavelength and fj{Xi) is 
the flux at time j and wavelength A^ . So as to better illus- 
trate the relative changes in each wavelength bin, we only 
plot the normalised standard deviation, i.e., a(Xi)/ f{Xi). 
This measure has the advantage of simplicity and univer- 
sality, but has the disadvantage that it makes disentan- 
gling facular and spot variability difficult. Secondly, we 
look at time series in a number of selected wavelength bins. 
These include the VIRGO/SPM filter bands, which allow 
us to compare our model, SORCE/SIM and VIRGO/SPM 
with each other, and a number of bands that stand out in 
the variability plots. 

Before calculating the final RMS spectra, the mean 
spectra and the time series, we removed obvious outliers 
in the SIM data. This was done for each wavelength bin 
individually, by removing data points that deviated from 



SIM detrended 
instrumental 




260 280 
wavelength [nm] 



Fig. 1. The normalised standard deviation between 220 
and 310 nm as derived from SIM/UV between 2004 April 
21 and August 1. The red and black lines show the nor- 
malised standard deviation of the smoothed data after 
the removal of outliers. The difference between the two 
variability spectra is due to the removal of a linear slope 
that reduces the variability for wavelengths above about 
260 nm as indicated by the black line (see section [ST] for 
a description). The dotted line traces the instrumental 
noise. Note that the instrumental noise exceeds the bi- 
nomially smoothed signal for wavelengths below approxi- 
mately 230 nm. 



the mean by more than ka. The cutoff factor, fc, was var- 
ied with wavelength to account for the different aspects 
of the faculae-dominated variability in the UV and the 
spot-dominated variability at longer wavelengths. We thus 
applied a symmetric cutoff in the UV, generally clipping 
data points more than 3.5a from the mean. In the visible 
and infrared, typically data points more than 2cr above 
and 4.5tT below the mean were clipped. The clipped data 
were replaced by median values of the two previous and 
subsequent exposures. 

The variability plots are shown in Figs [1] and O and 
will be discussed in the following sections. The plots also 
show the derived instrumental noise level. Note that SIM 
is generally not photon-noise limited, but analog-to-digital 
converter (ADC) limited with about 2 bits of noise on 
a 15-bit converter range. It requires a dynamic range of 
aproximately 100 to measure the signal, so for weak sig- 
nals the noise level becomes comparable to the solar vari- 
ability signature. Apart from these random noise contri- 
butions, additional residual systematic trends caused by 
the imperfect prism degradation and temperature correc- 
tions can still be present in the time series. As the analog- 
to-digital noise is essentially random, the application of 
a (non-phase shifting ) filter is appropriate. Here we use 
binomial smoothing ( Marchand &: Marmet 19831 ) in the 
time domain with two passes of the lowest-order (1,2,1) 
filter, meaning that we will be insensitive to variability on 
time scales shorter than about 1.5 days. 
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3.1.1. SIM/UV 

Solar variability is very much higher in the UV than in 
the visible and infrared, and robust measurements with 
variations of the ord er of several percent are expected. 
Harder et al. I (|2005bh have shown that the response of the 
UV instrument becomes more unreliable towards the blue 
end of the wavelength range. Fig. [1] shows the normalised 
standard deviation for the SIM/UV data. The red line 
indicates the normalised standard deviation of the original 
Version 10 data set, once outliers have been removed. The 
black dotted line indicates the instrumental noise. 

Over the time span considered here, the data from the 
UV detector show a slow, almost linear, decrease that in- 
troduces substantial variability and can be picked up in 
the 264 to 277 nm and 290 to 300 nm regions in par- 
ticular. This decrease could be either instrumental or in- 
strinsically solar in which case it would imply a slow de- 
crease of the normalised solar UV irradiance at the 5000 
ppm level over a 3-month time span. A solar origin is 
supported b y a comparison t o the SORCE/SOLSTICE 
instrument (jSnow et al.l 120051 ) over the same time span. 
While the SOLSTICE trend deviates in the first couple of 
weeks, it generally agrees with the SIM measurements for 
the remainder of the time. The red line in Fig. [T] shows 
the normalised standard deviation when the linear trend 
is removed. In this case, the variability around 270 nm 
decreases by a factor of 2 in better agreement with what 
is seen in the models (see Sec. 14. 4p . 

The measured standard deviation agrees well between 
the smoothed and unsmoothed data for wavelengths larger 
than approximately 260 nm where it also exceeds the in- 
strumental noise by more than a factor of two. For wave- 
lengths below 240 nm the instrumental noise becomes 
comparable to the data variability. This indicates that the 
measured variability is largely instrumental on the roughly 
12-hourly timescales considered here, and can thus be 
significantly reduced by binomial smoothing as indicated 
by the red line on Fig. [TJ Indeed, for wavelengths below 
235 nm, the variability of the smoothed data falls below 
the instrumental noise. 



3.1.2. SIM/VISl 

Fig.[2]shows the variability measured with VISl for wave- 
lengths between 310 and 550 nm. The red and black lines 
show the smoothed and unsmoothed data, respectively. 
Also shown is the instrumental noise level (dotted line). 
The figure illustrates that the instrumental variability in- 
creases dramatically for wavelengths shortward of approx- 
imately 400 nm (see also figure 2 in IWoodsll2007l ). While 
there is no marked difference between the smoothed and 
unsmoothed data above 400 nm, indicating that we mea- 
sure a predominantly solar signal, the variability of the 
smoothed data is lower (by a factor of about 1.5) at shorter 
wavelengths and falls below the instrumental noise level. It 
is thus not straightforward to estimate the solar variabil- 
ity from the data available here, or indeed even estimate 



SORCE/SIM 
SIM smoothed 
instrumental 
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Fig. 2. Variability between 2004 April 21 and Aug 1 as 
recorded with the SIM/VISl instrument. The black line 
shows the normalised deviation of the Version 10 data af- 
ter the removal of outliers. The red line is for the variabil- 
ity of binomially smoothed data, while the dotted line in- 
dicates the instrumental noise. Also shown is the variabil- 
ity measured in the VIRGO blue and green filters (green 
triangles). The short data stretch below 308 nm is as mea- 
sured with SIM/UV. 



unambiguously the range up to which the smoothed data 
represent solar rather than instrumental signal. 

Overall, an increase in the standard deviation is ex- 
pected for lower wavelengths, though the variability seen 
between 310 and 350 nm is clearly too high. We can take 
some guidance from the standard deviation of 500 ppm 
recorded with the SIM/UV detector at 300 nm. We would, 
however, caution against interpolating the variability be- 
tween 300 and 390 nm, despite the similar variability levels 
observed at both wavelengths. Not only does the region 
contain a number of intermediate-strength lines, it also 
coincides with the expected switch-over between the facu- 
lar and spot-dominated regime on rotational time scales. 
Depending on the exact balance between facular brighten- 
ing and sunspot darkening, both effects can almost cancel 
each other out. This would explain, e.g., why the variabil- 
ity is lower at 385 and 395 nm compared to 400 nm. 

Not shown on Fig [2] is the VISl variability above 
550 nm, as it is mainly featureless: it shows a slow de- 
crease between 550 nm to 800 nm where it ranges from 
350 ppm down to about 270 ppm. For longer wavelengths 
(> 820 nm) it shows an upturn. This var iability increase 
was already noted bv lHarder et al. I (|2005bh who attributed 
it to the incomplete removal of temperature-induced vari- 
ability in the instrument. The variability recorded by 
SIM/VISl is shown over the full wavelength range and 
discussed further in Sec. 14.41 where it is compared to the 
model results. 
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3.1.3. SIM/IR 



3.4. VIRGO 



SIM/IR, the infrared detector records the solar spectrum 
between 850 nm and 1.66 fim. The data in the IR suf- 
fer from occasional sudden data jumps in time. The data 
become particularly noisy at the detector edges. We thus 
only use data for wavelengths between 980 and 1600 nm 
in the following. In this wavelength range, the data are 
very uniform with a normalised standard deviation be- 
tween 230 and 300 ppm. 



3.2. SORCE/TIM 

With TIM, the Total Irradiance Monitor, the SORCE 
satellite also carries a solar radiometer to measure total so- 
lar irradi ance. The TIM inst rument has been described in 
detail bv IKopp fc Lawrence! lioOSi) a nd first results have 



been presented in iKopp et al. Ulooi). The instrumental 



noise level is less than 2 ppm and the instrumental stabil- 
ity is corrected to <10 ppm/yr, so the TIM data require no 
long-term gradient removal or high-frequency temporal fil- 
tering for the analyses here using Version 5 data. The TSI 
measured by SORCE/TIM is about 5 W m^^ lower than 
the TSI measured by ot her radiometers in sp ace, such as, 
ACRIM-III and VIRGO (i Frohhch et allll997[) . The irradi- 
ance changes of TIM, however, agree extremely well with 
those of the other radiometers, not only over the three 
months considered here, but also over the whole life time 
of the SORCE mission. Here we use TIM as representative 
for the TSI. As we consider relative changes in TSI only, 
we have normalised the modelled and SIM-integrated data 
(see SecO]) to the SORCE/TIM values. 



3.3. UARS/SUSIM 

The Solar Ultraviolet Spectral Irradiance Monitor 
(SUSIM) is a dual dispersion s pectrometer instrumen t 
that operated from 1991 to 2005 (jBrueckner et al.lll993l ). 
It was one of the 2 UV experiments on board UARS 
(Upper Atmosphere Research Satellite). SUSIM has been 
monitoring solar irradiance in the range from 115 to 
410 nm with a spectral resolution between 0.15 and 5 nm. 
We u se the daily level 3 BS V22 data with sampling of 
1 nm (jPlovd et al.1l2003H Floyd, priv . comm. 2006) avail- 
able at ftp: / /ftp.susim.nrl.navy.mil[ Calibration of the 
changing responsivity of SUSIM's working channel was 
done through a combination of measurements of four on- 
board deuterium calibration lamps and solar measure- 
ment s by less frequently exposed referen ce optical chan- 
nels (|Prinz et al.lll996l : IPlovd et al.lll998h . The long-term 
uncertainty of irradiance measurements (Icr) is about 2- 
3% at A > 170 nm, « 5% at A = 140 - 170 nm 
and increases to around 10-20% at s horter wavelengths 
( Woods et aDll996l: iFlovd et al.lll998ll2003bl) . 



The VIRGO/SPM instrument onboard SOHO measures 
solar variability in three wavelength bands, centred on 
402 (blue), 500 (green) and 862 nm (red) with band- 
widths (PWHM) of 5.4, 5.0 and 5.7 nm, respectively. 
The data presented here are level 1.7 daily averages and 
have been obtained from the SOHO data archive. SPM 
measurements suffer from strong and non-linear degrada- 
tion, so that stretches longer than one month need to be 
corrected carefully before they can be used for compari- 
son purposes. To correct for the degradation, we divided 
the VIRGO/SPM data by VIRGO TSI data and fitted 
a quadratic function to 9 data points that coincide with 
times of low solar activity. This essentially pins the long- 
term behaviour in the colour channels to that of the TSI 
during quiet-Sun phases. 

Note that a newer SPM data release has recently be- 
come available where most of the long-term degradation 
has been removed (Frohlich 2007, priv. comm). A compar- 
ison between our corrected data with the new data release 
shows that the variability amplitudes agree very well. A 
small amount of (possibly spurious) long-term variability, 
however, remains even in the new data set. Rather than 
carrying out a similar procedure as outlined above, we 
decided to use the old, but corrected data set. 



4. Comparisons of the SATIRE model to 
SORCE/SIM and SOHO/VIRGO 
measurements 

4.1. Total solar irradiance 

As a first test, we compare the modelled total so- 
lar irradiance to the SORCE/TIM measurements as 
well as to the SORCE/SIM 'total' solar irradiance. 
This SIM pseudo-TSI was obtained by integrating the 
(smoothed) SORCE/SIM measurements over the avail- 
able SIM wavelength range. As this does not cover the 
full solar spectrum, the resulting integrated irradiance is, 
at 1230 Wm-2, about 10% lower than the TIM measure- 
ments. This value is in reasonably good agreement with 
model expectations: we find that the model irradiance be- 
tween 220 and 1660 nm is 1207 Wm~^. The comparison 
between the time dependence of the modelled, SIM inte- 
grated and TIM measured irradiance is shown in Fig. [31 
Both, the SIM wavelength-integrated data and the mod- 
elled TSI were renormalised to match the absolute value 
of the TIM TSI. For the SIM-integrated TSI, renormali- 
sation should yield an upper limit for the variability am- 
plitude, as the missing part of the spectrum is mainly in 
the IR where variability levels are expected to be lower. 
We therefore also tried an approach whereby we added a 
constant offset. The true behaviour is then expected to 
lie between these extremes. We found both lightcurves to 
be very similar with no significant changes for the correla- 
tion coefficients, and therefore only present the normalised 
lightcurves in the following. 
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Fig. 3. Total solar irradiance (TSI) from May to July 
2004. The black diamonds linked by the solid black lines 
show the integrated SORCE/SIM data after binomial 
smoothing; the red triangles show the SORCE/TIM to- 
tal solar irradiance and the blue plus signs linked by the 
dotted lines indicate the integrated modelled irradiance. 
The model and SORCE/SIM data have been integrated 
between 220 and 1660 nm and have been normalised so 
that their mean matches the absolute value of the mean 
SORCE/TIM TSI. 

Table 1. A list of correlation coefficients for the differ- 
ent TSI determinations. The first two columns give the 
data sets used for the correlations, while the third col- 
umn lists the Pearson linear correlation coefficient r and 
its square. The fourth and fifth columns give the more ro- 
bust Spearman rank correlation, p, and its corresponding 
probability of a chance correlation. 



set 1 


set 2 


r [r^] 


P 


prob 




integrated SIM 


TIM 


0.97 [0.94] 


0.86 


9 X 10" 


26 


model 


TIM 


0.97 [0.94] 


0.89 


9 X 10" 


28 


integrated SIM 


model 


0.92 [0.84] 


0.72 


1 X 10" 


20 



The agreement between the SORCE/TIM, integrated 
SORCE/SIM measurements and the modelled TSI is good 
overall, as borne out by the correlatio n coefficients and 
chance probabilities listed in Tab. [1] (see Press et al. 19861 
for more detail). The difference between the model and 
the data is of the same order as the difference between 
the two data sets. Closer inspection of Fig. [3l however, 
shows that some inconsistencies remain. The model ap- 
pears least reliable between June 21 and June 30. This 
is mainly due to a lack of MDI magnetograms and con- 
tinuum images, and the poorer quality of those magne- 
tograms and images that are available. Furthermore, the 
model seems to overestimate the facular brightening asso- 
ciated with the spot passage in July. The reason for such 
an excess brightening could either be that our facular con- 
trast calculations are too high for large active regions, or 
it could arise through errors in the feature identification. 



e.g., if some of the spot/pore magnetic flux were wrongly 
attributed to faculae. This may occur in particular when 
active regions are near the limb. Finally, uncertainties in 
the LOS magnetic field correction can lead to the mag- 
netic field strength and hence the contrast of the faculae 
being overestimated. 

The main difference between the integrated SIM data 
and the other two data sets occurs during the passage of 
the two small spot groups in May and in the period just 
after the sunspot passage in June. Compared to either 
TIM or SATIRE, the integrated SIM data show a larger 
fiux increase before the passage of the first spot group, 
followed by a stronger flux decrease during the passage of 
the second spot group. It is not clear what might have 
led to this difference, as the data do not appear partic- 
ularly noisy or discontinuous. At the end of June, SIM 
fails to pick up the facular brightening after the sunspot 
passage. In fact, the whole period between the two SIM 
data gaps (around June 24 and July 10) shows a different 
behaviour than expected from the model or the TIM data. 
In Sec. 14. 2i we show that not all wavelengths are equally 
affected by this problem. Integration over a bluer wave- 
length stretch, e.g., one that excludes wavelengths above 
800 nm for VISl, produces a flatter response during that 
time. The lightcurve is otherwise very similar and the re- 
sulting increase in the correlation coefficient is very slight, 
so that it has not been plotted here. 



4.2. Comparisons between SIM, VIRGO and the 
SATIRE model 

Comparisons of VIRGO short-term spectral var i ations 
and o ur model have been prese nted bv lFligge et al. ( 1998 . 
2000l ) and lKrivova et al.l (|2003h . Here we extend this work 
and compare model calculations to SORCE/SIM as well as 
to the VIRGO/SPM irradiances. The VIRGO/SPM filters 
are narrow, to the extent that the FWHM of the green and 
red filters lie below the corresponding SIM spectral reso- 
lution. The use of a detailed filter profile is meaningless in 
such a situation and simply employing a single SIM wave- 
length channel leads to overly noisy narrow-band fluxes. 
In order to be able to include a larger number of wave- 
length bands and achieve better signal-to-noise ratios, we 
used wider, rectangular fllters, ranging from 490 to 510 nm 
in the green, and from 830 to 900 nm in the red. In the 
blue, we used the VIRGO/SPM bandwidth, as the num- 
ber of wavelength points covered by the blue filter is sig- 
nificantly larger than in the green and red; furthermore, 
extending the blue filter is difficult as there is not much 
clean continuum cither side of it. The filter widths and 
central wavelengths as applied to the different data sets 
are listed in Tab. [H 

As discussed in Sec. [3T2l the SIM/VISl detectors suf- 
fer from an imperfectly corrected temperature response 
above approximately 820 nm. This proved particularly 
troublesome in May 2004 and also during the data gap 
in early July. In order to exclude systematic effects aris- 
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Table 2. Central wavelengths and filter widths. For the 
VIRGO filters largest and smallest wavelength indicate 
the extent of the FWHM. The green and red filters for 
the model and SIM data are simple rectangular filters. 



filter 




VIRGO 




SIM/SATIRE 




Ac 


ran] 


je 


Xc 


range 


blue 


402.6 


400.0 - 


405.3 


402.6 


400.0 - 405.3 


green 


500.9 


498.5 - 


503.3 


500.0 


490.0 - 510.0 


red 


863.3 


860.5 - 


866.0 


865.0 


830.0 - 900.0 


reds 








775.0 


750.0 - 800.0 



ing through this, we also tested an alternative red channel 
(reds, for short red) that was obtained by integrating be- 
tween 750 and 800 nm. This reds channel suffers less from 
temperature-induced variability. Note, however, that the 
response in the two red channels is different. Our model 
calculations suggest that the variability in the shorter 
channel (reds) should exceed that of the original VIRGO 
channel by approximately 10%. Comparing the SIM data 
also suggests a larger response of the shifted reds channel 
with respect to the original VIRGO channel (by approxi- 
mately 6%), though this measurement is uncertain as, on 
account of the superimposed spurious variability, it has to 
be determined from a shorter data train. 

Comparisons between the VIRGO/SPM data, SIM 
and our model are shown in Fig. (4] We find that a number 
of days stand out in all filters. As already found for the 
TSI, the strongest discrepancy regarding SIM data occurs 
around June 25, just after the June sunspot passage. SIM 
apparently fails to pick up the facular brightening as the 
active region is near the limb; this is particularly salient in 
the blue and green filters. Note that the red filters mirror 
some of the problems seen in Sec. 14.11 namely an enhance- 
ment in mid May just before the first spot group appears 
and a rise at the end of June before the data gap. These are 
even more pronounced in the original red filter (data not 
shown here) and we thus conclude that they are largely 
due to an incomplete removal of the instrumental temper- 
ature changes that affect the longer wavelength regions 
particularly strongly. 

The largest difference between the model and both 
SIM and VIRGO data is the much larger facular bright- 
ening before and after the July sunspot passage. This is 
very noticeable in all three filters, as indeed also for the 
TSI (see Sec. 14. 1|) . By contrast, we find that the sunspot 
darkening measured by VIRGO and SIM agrees very well 
with the models, in particular in the green and blue filters; 
in the red band there is a slight tendency for the darken- 
ing to be overestimated. An exception to the good fit is 
the very small spot at the end of May where the model 
underestimates the sunspot darkening in all filters. 

One way to judge the tightness of the correlation be- 
tween two data sets is to consider histograms of their 
fractional differences as shown in Fig. [5l Because of the 
response problems of SIM/VISl for wavelength in excess 
of about 820 nm, we have used the short red filter (reds) 



to obtain the fiuxes for the SIM-to- model comparisons. 
The original central wavelength was used for the model-to- 
VIRGO comparison and the two different filters are used 
when comparing VIRGO to SIM. While the number of 
data points considered here is relatively small, we find that 
most of the histograms resemble Gaussians, some of them 
with noticeable skew. The largest deviations are seen for 
the blue filter for the SIM vs SATIRE comparison, where 
the distribution is very broad and may be double-peaked. 
In most of the cases, however, we can use the width of 
the standard deviation of a Gaussian fit to the histograms 
to characterise the scatter between the different data sets. 
These are listed in Tab. \3\ 

In order to quantify the fits further, we list the correla- 
tion coefficients in Tab. [31 To calculate the correlations, we 
binned the SORCE/SIM and model data onto the same 
grid as the daily VIRGO/SPM resuhs. Fig. [6] shows the 
correlation plots of the VIRGO/SPM and the SIM filter 
observations with respect to the model calculations. The 
red and black lines show the best-fit lines, assuming that 
all data sets suffer from equal relative errors. These errors 
were estimated to be of the order of 100 ppm in the green 
and red filter, and of the order of 180 ppm in the blue 
filter (see also Tab. [3]). The dashed blue line indicates a 
slope of unity as would be expected for a perfect match. 

When considering all filters together, the correlation 
coefficients and histogram widths suggest that the best 
agreement is found between the VIRGO data and SATIRE 
calculations, while comparisons fare least well for SIM ver- 
sus SATIRE. The correlation coefficients and plots indi- 
cate that the agreement between the model and the obser- 
vations is best for the green filter, but slightly less good for 
both the red and blue filters. In the following, we discuss 
the fits for the individual filters in more detail. 

The correlation coefficients in the red filter show a rel- 
atively large variation, ranging from a tight correlation 
with coefficient 0.96 (r^^O.OO) between VIRGO data and 
model calculations, down to a coefficient of 0.77 (r^=0.59) 
for SIM/red to model comparisons. This latter value can 
largely be explained by the temperature-induced sensitiv- 
ity problems for SIM/VISl above about 850 nm. A com- 
parison between the model and SIM data for the slightly 
shorter reds filter where this is less of an issue gives a sig- 
nificantly higher correlation coefficient of 0.92 (r^=0.84). 
Despite the filter shift and the associated change in re- 
sponsivity to spot and facular passages, there is also an 
increase in the correlation coefficient (from 0.89 to 0.96) 
when comparing VIRGO red with SIM/reds. 

In terms of the correlation coefficient, there is no sig- 
nificant difference between using the model calculations 
for the red or reds filter to compare with a given observed 
time series. This is expected, as the basic features of the 
model differ only very slightly between the two wavelength 
bands. The amplitude of the variability, however, matches 
much better when the model calculations are carried out 
for the filter appropriate for the comparison data. 

The slopes derived for the red correlations show that 
SATIRE overestimates the amplitude of the variability by 
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about 5% (SIM) and 10% (VIRGO). For the SIM analy- 
sis, we used the alternative 'short' red filter as this is es- 
sentially unaffected by the temperature effects. Note that 
the very small gradient of 0.88 derived for the VIRGO 
vs model comparisons is in part due to a single outly- 
ing model point on June 21, shown as the red triangle at 
(0.99935, 0.99975) in the bottom plot of Fig. M (see also 
Fig. H]). If this outlier is excluded, the gradient increases 
to 0.91 which is in better agreement with the results for 
SIM. A further reason for the low gradient could be that 
the red filter band includes the redmost line of the Ca il 
IR triplet that might not be modelled well by SATIRE. 

The correlations between all the data sets are excel- 
lent in the green filter. In fact, they are better than those 
determined for the TSI comparisons, presumably because 
they are limited to the similar narrow spectral region and 
therefore sample a well-definied region of the solar atmo- 
sphere. We also find that the gradients of the best-fit lines 
are near unity, indicating that the amplitude of the vari- 
ability agrees between all three datasets. In the blue fil- 
ter, the agreement between the model calculations and 
VIRGO is comparable to that of the TSI comparisons, 
though the model fares less well with respect to the SIM 
data. As indicated by the correlation gradients, the model 
appears to underestimate the variability by between 10 
and 20%. We consider the 20% derived from SIM data to 
be less reliable, mainly because the data follow a slight 
degradation-like long-term trend. This is again most ob- 
vious at the beginning of May. A correlation analysis with 
data after May 11 yields a slope of 1.11 with respect to 
the model, and agrees well with our findings for the blue 
VIRGO fiher. 



Table 3. Table listing the standard deviation to gaussian 
fits to the histograms (ci, in ppm) the linear correlation 
coefficient (rj), its square, and the regression slope (mi) in 
the three VIRGO/SPM filters. Golumns 2 and 3 list the 
quantities for the comparison for VIRGO/SPM against 
SORGE/SIM and the SATIRE models, respectively, while 
column 4 compares the SORCE/SIM measurements to 
the SATIRE model calculations. Standard deviations in 
curved brackets indicate that the histogram had signifi- 
cant skew and/or sidelobes and that a Gaussian was not a 
good fit. In the rows with the regression slopes, the square 
brackets give the errors on the regression slopes, assum- 
ing equal errors for SIM, VIRGO/SPM and SATIRE. Note 
that there are no VIRGO data for the redg filter; the last 
three lines in columns 2 and 3 thus give correlation coef- 
ficients and slopes comparing the VIRGO red band and 
the SIM or model reds bands. 







VIR vs SOR 


VIR vs SAT 


SOR vs SAT 


N 




84 


79 


71 


0"bluc 




172 


146 


(231) 


rbluc. 


'"bluJ 


0.95 [0.90] 


0.96 [0.92] 


0.89 [0.79] 






1.07 [0.04] 


1.10 [0.05] 


1.20 [0.06] 


^grccn 




109 


95 


92 


^grccn , 


[^grccn] 


0.98 [0.96] 


0.97 [0.95] 


0.96 [0.92] 


^^grccn 




0.78 [0.05] 


1.04 [0.03] 


0.99 [0.04] 


0"rcd 




(85) 


71 


(177) 


Trod, [r 


2 1 
rcdj 


0.89 [0.78] 


0.95 [0.90] 


0.77 [0.59] 






1.46 [0.04] 


1.14 [0.06] 




^rcd.s 




93 




91 


^red,s , 




0.96 [0.92] 


0.96 [0.91] 
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1.04 [0.06] 



4.3. Comparisons between SIM, SUSIM and the 
SATIRE model 

During 2004, the solar UV spectrum and its variability 
was also recorded with the UARS/SUSIM instrument. In 
this section, we compare these measurements with the 
SORCE/SIM measurements and the model calculations 
for wavelengths between 170 and 320 nm. Fig. [7] shows 
a plot of the normalised standard deviation for data and 
model calculations between May 1 and July 31 in 2004. In 
order to reduce confusion in the plot, we have binned the 
SORCE/SIM and UARS/SUSIM data in the wavelength 
domain before carrying out the variability analysis. The 
binning factors are detailed in the caption of Fig. \7\ A 
number of striking features are apparent in the plot, and 
are discussed below. 

The SORGE/SIM and SATIRE data show a large 
increase in variability at about 205 nm. This is most 
likely due to Al and Ga opacity edges between 200 and 
210 nm. A significant decrease of the solar brightness tem- 
perature around thi s wavelength was already observed by 
Widing et al. (1973) based on data from rocket flights. It is 
very noticeable, however, that the variability recorded by 
UARS/SUSIM is much lower than both the SORGE/SIM 



and the SATIRE variability While the UARS/SUSIM in- 
crease might appear weakened because of the higher (in- 
strumental) variability seen above 200 nm and the lower 
velocity resolution, we expect the SUSIM data to best re- 
flect the solar variability below 200 nm during the time pe- 
riod considered here. The main reason for the (excessive) 
increase in variability of the SATIRE model is the break- 
down of the LTE assumptions and the use of opacity distri- 
bution functions (ODFs) rather than detailed line-opacity 
calculations. In the case of SIM, the jump in variability 
is not surprising either, as its detector is not expected to 
perform well at these wavelengths (see Sec. 13.1. T|) . Better 
results should be provided by SORCE/SOLSTICE. 

In the wavelength region between 210 and 290 nm, the 
agreement between the model calculations and the SIM 
measurements is varied. While some features such as the 
Mg II h&k lines (280 nm), and the regions from 220 to 
232 nm, from 255 to 270 and above 290 nm match well, 
other wavelength regions show large disagreements, see, 
e.g., the lines of Mg i at 285 nm, or the complex sets of 
lines around 240 and 250 nm. As above, these disagree- 
ments are due mainly to the assumption of LTE and the 
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use of ODFs (jKuruca 119921 ). Uncertainties in the model 
atmospheres also contribute, though probably to a much 
smaller extent, since the differences are largest at the 
wavelengths of strong lines showing strong NLTE effects. 
Thus the difference in the behaviour of the Mg i and Mg ii 
resonance lines can be explai ned quite well if NLTE effects 
are taken into account fUit enbroek fc Briandl 1995 ). 

Overall, the agreement between SIM and SUSIM is 
reasonably good between about 210 and 290 nm. While 
the variability recorded with SUSIM is higher due to its 
lower sensitivity, the features recovered agree well and the 
measured relative variabilities are not too different. Above 
290 nm UARS/SUSIM shows a relatively poor response 
that swamps the solar variability on mid to short-term 
time scales. The responsivity and noise characteristics of 
SUSIM have bee n well documented and are discussed, e.g., 
(|l99d ). 



4.4. Comparison of SIM with the SATIRE model over 
the whole wavelength range 

The main advantage of SIM over the VIRGO/SPM chan- 
nels is that a much more complete sampling of wavelengths 
is available. Fig.[H]shows the RMS variability between May 
and August 2004 over the whole analysed SORCE/SIM 
wavelength range. The normalised standard deviation of 
the SORCE/SIM data is indicated by the solid black lines 
in panels a to c, though note that the grey-shaded area on 
panel a indicates the range in variability that is obtained 
depending on whether a linear trend is removed from the 
data or not. The horizontal bars in Figs. [8] b and c in- 
dicate the wavelength regions where the measured vari- 
ability is dominated by instrumental noise (black dotted 
lines). The modelled variability is indicated by the blue 
lines. The bottom plots (d, e and f) show the contribu- 
tions of the spots (dotted black lines) and faculae (dashed 
purple lines) to the overall variability. To calculate these, 
we replaced the facular (resp. sunspot) contribution by a 
quiet-Sun contribution. The curves show very clearly that 
the wavelength dependence of spot variability is spectrally 
much smoother than the facular variability. This has to do 
with the fact that the darkening due to spots is dominated 
by the drop in continuum intensity. Changes in spectral 
lines produced by the lower temperature in spot umbrae 
and penumbrae play a secondary role. For faculae, the 
absolute temperature difference is less pronounced (espe- 
cially in the lower atmosphere) , and it is the different tem- 
perature gradient that produces changes in the continuum 
as well as in the lines. Especially at shorter wavelengths, 
individ ual and groups of lines provid e the dominant cont ri- 
bution (|Mitchell fc Livingstoniri99ll : IUnruh et al.ll2000( ). 

Below about 280 nm, the variability due to faculae 
generally exceeds that due to spots; then follows a region 
up to 400 nm where they are mostly of comparable mag- 
nitude. Above 400 nm, the modelled spot variability is 
always larger than the facular variability, and the com- 
bined variability follows the spot variability closely. Note 



that there is a further cross-over around 1.4 /im where the 
modelled variability of the spots and faculae drops below 
the total variability. This marks the transition where the 
facular model becomes dark averaged over the solar disk 
and thus no longer acts to counterbalance the spots. The 
fact that the model shows a smaller variability than SIM 
at these wavelengths cannot be due to this property of the 
model; dark faculae enhance the darkening due to spots, 
and thus increase the standard deviation (see also Sec. [5]). 
In other words, if the model faculae were bright at these 
wavelengths the discrepancy between the SATIRE model 
and SIM data would be larger. 

The shift in importance away from faculae to spots 
at wavelengths as low as about 300 nm is expected when 
considering variability on the order of a couple of month, 
i.e., on the rotation time scale when the influence of spots 
tends to dominate the TSI. On longer time scales such as 
that of the solar cycle, however, bright small-scale mag- 
netic features dominate the TSI variations, and are thus 
also expected to dominate variability in both the UV and 
visible. 



4.5. SIM & SATIRE time series 

In this section, we compare the variability in a number of 
wavelength bands in more detail. The wavelength bands 
have been picked so as to show the change in behaviour 
going from the UV around 220 nm, up to the IR at 1.5 /xm. 
They are typically also chosen at wavelengths where the 
relative variability is high and the data quality is good. 
The wavelength bands are listed in Tab. |4l along with 
their band widths and the number of data points included 
in the integration of the SIM data. So as to improve the 
S/N level of the resulting time series, the UV bands, in 
particular, have been chosen to include a large number of 
wavelength points. Note that the noise level of the mod- 
elled time series is largely governed by the noise in the 
magnetograms and cannot be decreased by increasing the 
number of wavelength points considered; the absolute flux 
level, however, is influenced by the coarseness of the wave- 
length grid. 

The measured and modelled timeseries are shown in 
Figs [9] and [10] and the correlation coefficients as well as 
the slopes between the model and the data are listed in 
Tab. m The upper panels show the measured and mod- 
elled solar variability while the bottom panels show the 
modelled contributions of the faculae and spots to the ir- 
radiance variations. The figures illustrate again the rapid 
decrease in the variability towards longer wavelengths. 
Additionally, they show the change in the lightcurve as- 
pects as one moves from the UV to the visible. In the UV, 
as illustrated by the 230 and 280 nm bands (Fig. [9]), the 
influence of the spots is so small that their darkening effect 
is more than compensated for by the faculae, even on solar 
rotational time scales. Furthermore, the facular contrast 
increase towards the limb is not sufficient at these wave- 
lengths to counteract the projection effects. Consequently, 
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Fig. 8. Normalised standard deviation over the whole SORCE/SIM wavelength region. From left to right, the panels 
show the variability between 220 and 320 nm, between 290 and 600 nm and between 600 and 1600 nm. The top figures 
show the modelled and the observed variability in blue and black, respectively; in the left-most plot, the grey-shaded 
area indicates the range of variability for SIM, depending on whether a linear trend is removed. Also shown is the 
instrumental noise limit (dotted lines) and the variability measured in the three VIRGO filters (green triangles). 
The horizontal solid lines offset from the plots indicate the wavelength ranges where instrumental noise and artefacts 
dominate the variability. The bottom plots show the modelled variability, distinguishing between total variability (blue 
solid line), the spot (dotted black) and facular (dashed purple) contributions. 



the Sun appears brightest when the main spot groups are 
nearly at disk centre. 

At our resolution, this behaviour is no longer observed 
at longer wavelengths, such as in the Ca ii and the G 
bands. The spot darkening is now sufficient to offset the 
facular brightening, at least when the active regions are 
near disk centre. The overall lightcurves thus appear some- 
what confusing, with rapid sequences of peaks and dips, 
and very few stretches of 'quiet'-Sun behaviour. Further 
high-resolution calculations would be required to check 
whether this behaviour also holds for the individual line 
cores, or whether we are currently seeing a mix of line 
and continuum behaviour. The low wavelength resolution 
of the model and the failure to calculate exact line pro- 
files partly explains the relatively large difference in the 
observed and modelled mean flux in the Mg ii and, in par- 
ticular, 395 nm (Ca H&K) bands, where the model only 
includes 4 wavelength points. A further reason for the dif- 
ference in flux is due to our assumption of LTE, most 
markedly for Ca H&K. 

The situation again changes when looking at (contin- 
uum) bands and longer wavelengths in general (see left- 
hand panel in Fig. I10|) . The facular brightenings are now 
much weaker and mostly show a double-peaked aspect, 
i.e., the faculae produce most of the brightening when 



near the limb, and very little or even no brightening when 
at disk centre. Combined with the spot contribution, this 
leads to the familiar spot-dominated light curves, with 
small brightenings just before and after spot passages. 
Such behaviour is indeed seen for the TSI, and has been 
discussed at length in the previous section. 

In the NIR, finally, the facular brightenings become 
very weak, and might even disappear completely for the 
facular model atmosphere we employ. This is suggested by 
the model calculations for the 1550 nm band that differ 
markedly from those for the 1065 nm band. Fig. [10] sug- 
gests that the spots are not sufficiently dark at 1550 nm 
and indicates that the temperature of the spot model at- 
mosphere is too high in the deeper layers (1550 nm is close 
to the opacity minimum and thus carries information on 
the deepest observed layers). But while the model cal- 
culations appear to underestimate the flux decrease due 
to most of the active regions, they agree very well with 
the timing of the flux decrease; note that the dips due to 
the spot passages are significantly wider at 1.55 than at 
1.07 fim in the observed as well as in the modelled data. 
In the model, the wider spot passages are a consequence 
of the very small contrast of the faculae at that wave- 
length. Fig. [TT] shows a comparison of the model facular 
contrast for the six longer wavelength bands. The contrast 





Fig. 9. Top plots: flux variations for the 230 nm and the Mg II band (left-hand side) and of the 395 nm and the G 
band (right-hand side). The exact band widths are listed in Tab.[31 The diamonds linked by the solid lines represent 
SIM data, the plus signs and dotted lines show the model calculations. Note that the modelled timeseries have been 
binned onto daily values; the SIM data remain on their original time resolution, but have been smoothed using a 
binomial filter. Bottom plots: Modelled time series for the facular and spot contributions. The upper lines are for the 
faculae, the lower for the spots. The band wavelengths for the solid and dashed lines are indicated on the plots. 



at 1550 nm shows the lowest values throughout and be- 
comes negative near the disk centre (/x > 0.5). This low 
contrast means that only very little facular brightening 
is seen near the limb and thus leads to an earlier onset 
of the spot-induced darkening, as illustrated on the right- 
hand plots of Fig. \W\ 



We note that at 1.55 /im, there is rather poor agree- 
ment between the model calculations and the SIM data 
during May 2004. In particular, we find that the SIM data 
show a reversal in the relative spot strength during May. 
The first spot (centred around May 15th) is significantly 
stronger than the spot at the end of May in all wavebands, 
except at 1.55 iim, where the second spot appears darker. 
This could indicate that the temperature gradient in the 
two spots is different, although we cannot exclude uncor- 
rected data fluctuations. 



5. Discussion and Conclusions 

We have presented and compared SATIRE model calcu- 
lations and measurements of spectral solar variability on 
rotational time scales. The data and calculations cover a 
3-month time span from May to July 2004. In addition, 
we also compare modelled and observed time series of the 
total irradiance variability and the variability in a num- 
ber of selected wavelength bands. Such comparisons are 
particularly timely as SORCE/SIM is able to provide un- 
precedented observations over most of the range starting 
in the UV at approximately 220 nm and including the 
visible as well as the near infrared up to 1.6 /im. 

We find excellent agreement between the modelled 
total solar irradiance variations and the SORCE/TIM 
measurements. The absolute value of the wavelength- 
integrated SORCE/SIM measurements is in line with the 
expected model fluxes, and its variability agrees well ex- 
cept on a small number of days when the data quality 
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Fig. 10. Same as Fig. [9l though this time showing wavelength bands centred at 511 nm and at Ha respectively 
(left-hand side), and two near-IR bands centred at 1065 and 1550 nm, respectively (right-hand panels). 



was poorer (see Fig. [3]). We find correlation coefficients 
of 0.97 and 0.92 when comparing the modelled total solar 
irradiance with TIM and the wavelength integrated SIM 
measurements, respectively. 

The modelled and measured spectral variability over 
the three months is summarised in Fig. [8] for wavelengths 
between 220 and 1600 nm. Overall, we find good agree- 
ment between the model and the observations. Agreement 
is particularly good between 400 and 1300 nm. In the 
UV, where we also compare the SIM measurements to 
UARS/SUSIM, the agreement is somewhat patchy; some 
strong individual lines, such as the Mg ii h&k doublet, 
match very well, others, such as Mg I and Ca ii H&K, 
agree only poorly. This is not too surprising as we 
use opacity distribut i on an d assume LTE throughout. 
Uitenbroek fc Briandl (|l995l ) have shown that NLTE ef- 
fects can explain much of the different behaviour of the 
Mg I and Mg ii resonance lines. We also note that the 
resolution of our calculations is insufficient to resolve even 
the strong lines and to capture their complex behaviour. 
The role of spectral resolution in the context o f line v ari- 
ability has been discussed, e.g., in White et al. ( 2000f ). 



In the wavelength range between approximately 
310 nm and 350 nm, possibly even up to 390 nm, the 
response of both SORCE/SIM and UARS/SUSIM is too 
poor to determine solar variability on the rotational time 
scale. The best estimate of variability at those wavelengths 
is currently provided by the SATIRE model. The model 
calculations allow us to isolate the facular and spot con- 
tribution. This, together with the light-curves, illustrates 
very clearly the change from facular dominated variability 
at short wavelengths to spot-dominated variability above 
approximately 400 nm. 

In the visible, the observed and modelled irradiance 
variability matches well, though the decrease in the vari- 
ability at longer wavelengths appears somewhat steep in 
the model compared to the observations. We find, e.g., 
that the SATIRE model overestimates the variability be- 
tween about 600 and 800 nm by up to 20 % compared to 
the SIM measurements, while it underestimates the vari- 
ability around 1.5 /im by a similar amount. Note that for 
wavelengths between 800 and 1000 nm, the SIM detectors 
suffer from temperature-induced variability that cannot 
yet be fully compensated for; we were therefore unable to 
carry out meaningful comparisons at those wavelengths. 
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Table 4. Table listing the different bands used to calculate the time series. The first column gives the label as used 
for the plots, the second and third columns give the start and end wavelengths for the bands, while column four gives 
the number of flux points over which the integral was carried out for the SORCE/SIM data. Columns 5 and 6 give the 
mean fluxes (in units of W m~^) of the band as derived from the model calculations and the SIM measurements. The 
last three columns list the correlation coefficients, the gradient and the y-axis flux offsets of a linear fit of the data to 
the model; the gradient and offset for Ca ii is in brackets as the correlation coefficient is rather low. 



band 


As 


\e 


number of 


flux 


flux 


correlation 


gradient 


offset 




[nm] 


[nm] 


flux points 


(data) 


(model) 


coef 


(vs model) 




[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


230 


220 


240 


271 


1.02 


0.98 


0.90 


0.99 


0.06 


Mg II h&k 


277 


283 


38 


1.20 


0.71 


0.96 


1.03 


0.47 


Ca II H&K 


391 


398 


14 


7.21 


5.60 


0.51 


(0.95) 


(1.89) 


G-band 


420 


435 


22 


23.5 


21.3 


0.77 


1.35 


-5.3 


511 


507 


516 


7 


13.5 


15.3 


0.95 


0.88 


0.1 


Hq 


644 


668 


10 


32.0 


33.8 


0.96 


0.84 


3.6 


1065 


1050 


1080 


6 


16.7 


15.8 


0.89 


1.24 


-2.9 


1550 


1527 


1583 


10 


12.9 


14.6 


0.93 


1.33 


-6.4 



As a cross-check, we further compared the SATIRE and 
SORCE/SIM variability with the VIRGO/SPM measure- 
ments at 400, 500 and 860 nm. We found that the corre- 
lation coefficients for the model-to-data comparisons are 
typically very similar to those obtained for the SPM-to- 
SIM data comparisons. 

Our model suggests that the overall effect of faculac at 
1.6 ^m is one of darkening, though they appear bright at 
all wavelengths when seen close to the limb. Near 1.6 /im, 
the small brightness enhancement seen for faculae at the 
limb, however, is typically offset by the spots, so that 
active-region passages produce longer-lasting brightness 
dips at this wavelength than at shorter wavelengths. This 
is illustrated on th e right-hand panel of Fig. [TOl Contrary 
to Fontenla et al. 12004), we find no evidence for over- 
all bright factdae during the comparatively quiet period 
analysed here. An unambiguous contrast determination is 
difficult, however, as most large facular regions tend to be 
accompanied by dark spots whose exact contrasts are also 
unknown. 
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Fig. 4. Comparisons between detrended VIRGO/SPM 
data (red triangles and hnes), and SORCE/SIM (black di- 
amonds) as well as model data (blue plus signs linked by 
dotted lines) integrated according for the blue, green and 
red VIRGO filters. In the bottom plot, the SORCE/SIM 
data are for the redg filter (see text). Note the different 
scales for the y-axis on the three plots. 
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Fig. 5. Histograms illustrating the differences between the 
measurements and the model in the three VIRGO/SPM 
channels. The black and red histograms compare SIM, re- 
spectively SPM, against SATIRE. The blue dashed lines 
trace the histograms for SPM vs SIM. Note that bin 
widths for the blue filter have been doubled; and that the 
shorter reds filter was used for the SIM comparisons. 



Fig. 6. Plots of VIRGO/SPM and SORCE/SIM vs 
SATIRE in the blue (top), green (middle) and red (bot- 
tom) VIRGO channels. The black diamonds are for SIM 
vs model, the red triangles for VIRGO vs model data. The 
blue thin dashed line indicates a unit gradient, while the 
solid black and red lines show the gradients for the best 
linear fits for the SIM and VIRGO data, respectively. 
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Fig. 7. Normalised standard deviation for SORCE/SIM 
(black line), UARS/SUSIM (red) and the SATIRE model 
(blue) calculations. To lessen the confusion of the plot, the 
SIM data have been wavelength binned by factors of 10 
and 5 below and above 240 nm respectively; the SUSIM 
data were binned by a factor of two for wavelengths above 
210 nm. 
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Fig. 11. A plot of the modelled contrasts in the 6 longer 
wavelength bands from Tab. |31 The contrast decreases 
very strongly with wavelength, becoming negative at the 
disk centre at the near IR, as illustrated by the bottom 
curve for 1550 nm. The contrasts in the Mg ii and 230 nm 
bands are much higher and are not plotted here. Note also 
that these contrasts are maximum values that are scaled 
by the facular filling factors. 



